Pure-state quantum trajectories for general non-Markovian systems do not exist 
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Since the first derivation of non-Markovian stocliastic Sclirodinger equations, tlieir interpretation 
lias been contentious. In a recent Letter [Phys. Rev. Lett. 100, 080401 (2008)], Diosi claimed to 
prove that they generate "true single system trajectories [conditioned on] continuous measurement" . 
In this Letter we show that his proof is fundamentally flawed: the solution to his non-Markovian 
stochastic Schrodinger equation at any particular time can be interpreted as a conditioned state, 
but joining up these solutions as a trajectory creates a fiction. 
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It is well recognized that the continuous measurement 
of an open quantum system S with Markovian dynamics 
can be described by a stochastic Schrodinger equation 
(SSE). The pure-state solution to such an equation over 
some time interval, a "quantum trajectory" can be 
interpreted as the state of S evolving while its environ- 
ment is under continuous observation (monitoring). This 
fact is of great importance for designing and experimen- 
tally implementing feedback control on open quantum 
systems 0, 0, 3- If this interpretation could also be 
applied to non-Markovian SSEs 0], then this would 
be very significant for quantum technologies, especially 
in condensed matter environments, which are typically 
non-Markovian 

Previously we have argued that non-Markovian SSEs 
(NMSSEs) cannot be interpreted in this way The 
solution at any particular time can be interpreted as 
the system state conditioned upon some measurement of 
the environment [1]. But connecting up those solutions 
to make a trajectory is a fiction akin to trajectories in 
Bolimian mechanics Restricting to standard quan- 
tum mechanics, the basic problem is that for the state 
of S to remain pure, the bath field must be continuously 
observed to disentangle it from the system. For Marko- 
vian dynamics, this is not a problem, because the moving 
field interacts with S and, having interacted, moves on. 
But for non-Markovian dynamics, the field comes back 
and interacts again with S. Thus monitoring the field 
will feed disturbance back into the system, changing the 
average evolution of the state of S. That is contrary to 
the derivation of the NMSSE, which is constructed so as 
to reproduce on average the no-measurement evolution 
of S*. 

Recently, Diosi rcderived one form of NMSSE from a 
different starting point, and claimed that, contrary to 
the conclusions of Ref. Q, this allows an interpretation 
of the solutions as "true single system trajectories [con- 
ditioned on] continuous measurement" Q . Here we show 
by general argument, and an explicit calculation, that 
this claim is incorrect, and that the reformulation does 
not alter our earlier conclusion. 



The non-Markovian system. — Diosi considers a bath 
comprising an infinite sequence of von Neumann appa- 
ratuses An, each described by position and momentum 
operators Xn, Pn, n G {1,2, . . .oo}. (For clarity, we are 
using slightly different notation from Ref. Q.) The sys- 
tem interacts with the bath via the coupling Hamiltonian 



V 



S{t - Tn)Xpn, Tn 



en, 



(1) 



where X is an Hermitian system operator. Here the ex- 
plicit time-dependence plays the role of the free propa- 
gation of a bath field. This would seem to be a recipe for 
generating Markovian evolution, since S interacts only 
once with each An, which thus plays a role analogous to 
a small segment of a Markovian bath field. The novelty of 
Diosi's approach is to generate non-Markovian evolution 
by having the {Ak}'^^i prepared in an entangled state 
\(t>o)- In the position representation it is given by 



({2;fc}^il0o) oc exp 



-e'^y^^xiXmajri 

1.771 



(2) 



The continuum-time limit is e ^ 0, where the system 
is subjected to infinitely frequent, but infinitcsimally 
strong, interactions with the apparatuses. In this limit, 
a{t) plays the role of the correlation function for the 
bath. It is a real and symmetric function [12], and 
equals g^S{t) in the Markovian case. Assuming the sys- 
tem is initially in a pure state also, the Hamiltonian p]) 
produces an entangled system-bath state ]^'(r+)) imme- 
diately after the nth interaction. 

Diosi first considers the case where, immediately af- 
ter each time r„, the observable Xn is measured, yielding 
result Xn- This gives an unnormalized state for the con- 
ditioned quantum system, p(r^; {a^/ULi); given by 

Tr{A™}~.„,, [({^;}r=ii*K+))(*K+)i{^;}r=i)] . (3) 

with Tr [y5(r^; {x/}"^j^)] being the probability for the 
record {a^;}JLi- In the limit e ^ 0, this state (if ap- 
propriately scaled) will have a continuous but stochastic 
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evolution through time. The measurement of observable 
Xn does not disturb the future evolution of S because An 
never interacts with S again. Thus, there is no difficulty 
with interpreting this stochastic evolution as the trajec- 
tory of an individual system, with the average state at 
time t 



Pit) 



dxo ■ 



dx„ (4) 



being identical with that obtained simply by tracing over 
the bath (the apparatuses). 



(5) 



It is obvious, however, that p{t', {xi}il/^^) is not the 
solution of a SSE, for the simple reason that the state 
is mixed, not pure, even if it begins pure The 
mixedness arises because the interaction of S with An 
entangles S with Am for m > n, because initially An 
and Ajn are entangled. That is, the system becomes en- 
tangled with apparatuses that arc not yet measured. A 
mixed conditional equation state is not unexpected for 
non-Markovian systems. It has previously been shown in 
Rcfs. and [lH that it is possible to derive a mixed 
state quantum trajectory equation that reproduces the 
non-Markovian evolution on average by adding to 5 a 
fictitious system F, with the latter coupled to a moni- 
tored (Markovian) bath. A mixed state for S arises when 
the partial trace over F is performed. See Ref. [l^l for a 
comparison of this method with that of the NMSSE. 

The non-Markovian SSE and its interpretation. — The 
only way to obtain a pure state for S at time t is by 
measuring all the apparatuses with which the system is 
entangled. Specifically, Diosi shows that it is necessary 
to measure the set of bath observables {^s) : s <E [0, i]}, 
where z{s) is the "retarded observable" 



oc 

z(s) = 2e^Xka{s - Tk). 

fe=i 



(6) 



This is of course a different observable at different 
times s. The state conditioned on the result Zt = 
{z(s) : s G [0, t]} of this measurement at time i is a func- 
tional of z{s) for < s < i, which we will write as 
|V!t[Zt]). Diosi shows that this state is pure, and that 
it is the solution of the NMSSE 

d\MZt]) 



dt 



Xt[z{t)~2 I a{t^s)j-^ds] \ijt\Zt]) 



(7) 

Here, Diosi is working in the interaction picture with 
respect to the system Hamiltonian H] hence, the time 



dependence of Xt = e*^*Xe Equation ([7]) was first 



derived in Refs. [a, , but is very similar to that derived 
earlier in Refs. [5|. The ensemble average of solutions of 
this NMSSE reproduces the reduced state of the system: 



Here in taking the expectation value, z{t) must be treated 
as a Gaussian noise process with correlation function 
E[z(t)z(s)] ~ a{t — s), as appropriate for |^o)- This con- 
vention is indicated by the notation ijj (as oppiosed to ip) 
for the state. 

The contentious issue is not whether the solution 
has an interpretation in standard quantum me- 
chanics. As just explained, this state is the conditioned 
state of S at time t if an all-at-one measurement of the 
set of bath observables {z{s) : s G [0,t]} were made at 
that time, yielding the result Zt- The contentious issue 
is: can the family of states |^t[-2^t]) for < t < 00 be 
interpreted as a trajectory for the state of a single sys- 
tem, conditioned on monitoring of its bath. Diosi claims 
that it can be so interpreted, and that the required mon- 
itoring is simply to measure z{to) at time Tq , z{ti) at 
time Ti and so on. At first sight this monitoring may 
seem equivalent to the all-at-once measurement described 
above. But in fact it is not, as we will now explain. 

A measurement of z(t) at time t~^ involves measuring 
apparatuses that have not yet interacted with S. This is 
necessarily so because the symmetry of q;(t) means that 
z{t) contains contributions from Xm for some Tm > t 
(except for the Markovian case of course). Consequently, 
z{t) does not commute with pm for some > i, and 
the measurement will therefore disturb these momentum 
observables. But these are precisely the observables that 
will couple to the system via ((U, and thereby disturb it. 
Thus, as soon as the first measurement is performed, of 
z{to) at time tq, S ceases to obey the NMSSE. Whatever 
stochastic evolution it does undergo, it will not reproduce 
the reduced state of the unmeasured system pit). 

It might be thought that it would be possible to avoid 
this alteration of the future evolution of the system by 
repreparing the apparatuses Am for Tm > t in their pre- 
mcasurement states. However, this is not possible; before 
the measurement, these Am were entangled with the sys- 
tem S and the other apparatuses. The correlation of 
these Am with S is why the system state p(r„; {xi}^^^^), 
conditioned on measuring the apparatuses after they 
have interacted with the system, is mixed. The evolu- 
tion of this state over time is the only true quantum 
trajectory for a single system, and its mixedness is an 
inevitable consequence of the non-Markovian dynamics. 
In fact, we now show by explicit calculation that the 
monitoring Diosi suggests does not even produce pure 
conditioned states of S' — it also leads to mixed states. 

A simple example. — We consider the case where the 
bath consists of two apparatuses and e = 1 . Thus there 
are just three relevant times, tq = (the initial time), 
Tj^ = 1 (just after the interaction with Ai) and T2 ~ 
2 (just after the interaction with A2). Without loss of 
generality, we can write the initial Gaussian entangled 
state of the bath, analogous to Eq. as 



p{t)^E[\MZt]){MZt 



(8) 



(t)o{xi,X2) = cexp[— (x^ -t- ^2 + 2axiX2)], 



(9) 
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where = 2\J\ — a? /-k. Here < a < 1 parametrizes 
the initial entanglent between the apparatuses. The ana- 
logue of Eq. ^ defines two operators, 



zi — 2{xi + 0x2), Z2 = 2{x2 + ail). 



(10) 



Let us consider the unconditioned evolution of the sys- 
tem. At the initial time tq the total state is 

l^'o) ^ J^4'o{xi,X2) \xi)^ 1x2)2 \ipo)dxidx2, (11) 

where the final ket (with no subscript) denotes a state of 
S, and the subscript on the integral sign indicates it is a 
double integral. This evolves to the following state im- 
mediately after the interaction with the first apparatus: 



1*1) 



(l)oixi,X2) \xi 1X2) 



X \Xi) {Xi |-0o) dxidx2dXi 



(12) 



Here \Xi) denote eigenstates of Xi = X{ti), which we 
have taken to have a continuous spectrum for simplicity. 
Finally, after the second interaction, the total state is 

1*2) = J MXUX2)\X1+XI)^\X2+X2)2\X2) 

X {X2\Xi){Xi dxidx2dXidX2. (13) 

From Eq. p^ . the reduced state for the system at time 
T2 is simply 



P2 



Xi — X[ X2 — X2 



1X2) {X2\Xi){Xi\ 



4 \ 2 ' 2 
X {i:\X[) {X[\X!2) {X!2\dXidX2dX[dX!2. (14) 



All-at-once measurement at time r^. 



It is convenient 

to use, rather than the observables z„ pO|) . the scaled 
observables 



yi = zi/2 = xi+ ax2 = Ci(ii,i2), 
m = ^2/2 = X2 + axi = C2(ii, £2). 



(15) 



A measurement of z„, or y„, is described by the projector- 
density n„(i/„), defined by 

n„(y„) = J dxi J dx2 7ri(a;i)®7r2(a;2)^(j/n-CK(a;i,a;2)), 

„ (16) 

where 7r„(x) = \x)^ Note that, unlike 7r„(x), Tln{y) 
is not a rank-one projector; it is in fact a rank-infinity 
projector. It satisfies / dylln{y) = 1 and II„(?/)II„(y') = 
S{y—y')Iln{y) (no sum over n implied). It is obvious from 
the definition (|10|) that the two measurements commute. 

Consider first the case where at time T2 projective 
measurements of yi and y2 are performed. This yields 



*2 (2/1,^2; 



n2(y2)ni(yi)|*2) 



Vi - ay2 



1 -a2 



2/2 



1 - a 



^2(yi,y2) 



(17) 



where the conditional system state "02(2/1, 2/2) ) is 



c y^exph(Xi-yi)^-(X2- 2/2)1 

a^(2/? +2/1) - 2a2/i2/2 



X cxp y—2aXiX2 

X {X2\Xi){Xi\4,o)dXidX2 



1 -a2 



1X2 



(18) 



Obviously 5* is no longer entangled with {^1,^2}- This 
is as expected since the operators yi and y2 are linearly 
independent, and jointly measuring these is equivalent to 
jointly measuring xi and X2- That is, the measurement 
at time T2 effects a rank-one projective measurement on 
the bath, disentangling it from the system. Moreover, it 
is easy to verify that, as expected. 



1 



■*/'2(2/i,2/2)) (■(/'2(2/i,2/2) dyidy2 = p2- (19) 



This establishes that Eq. p8)) is indeed the discrete-time 
analogue of the solution of the NMSSE ([7]) at the relevant 
time (here r^). 

Monitoring (measurements at and T2 )■ — Now con- 
sider the case that Diosi claims is equivalent to the above, 
namely measuring yi at time and 1/2 at r^. From 
Eq. (fT2|) . the conditional total state at time is 



*i(2/i)) = ni(2/i)|*i) 
e-(i-')^ 



(20) 
V'i(2/i)) , 



I2/1 — ax)^ \x)2 dx 



where the conditional system state is 

Myi)) = cexp [-(yi - Xif] |i^o) ■ (21) 



So far we have a pure state for the system, as expected 
from Diosi's argument. However, at the very next step it 
breaks down. Because the measurement of the bath has 
disturbed it, we cannot use the state ((T3|) to calculate 
the next conditioned state. Rather, we must calculate 
the effect of the interaction between S and A2 on state 
The new entangled system-bath state at is 



(1-')-' \y^-ax)^\x + X2)2dx 
X 1X2) (^2 Myi))dX2 (22) 



*2|l(2/lj/ = /^^ 



Here the 2|1 subscript indicates that the state is at time 
T2 but the measurements it is conditioned upon was per- 
formed at time . 

After the second measurement we have 



^211,2(2/1,2/2)) = 62(2/2) ^'211(2/1) 



(23) 
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which evaluates to 

+ yi- ay2 



112 - a?X2- aiji 



X exp 



{X2 + ayi - 2/2)' 



l-a2 

X 1X2) (X2I exp[-(Xi - yi)^] IV-o) dX2 (24) 

Note that this is an entangled state between S and the 
bath — it is not possible to define a pure conditional state 
for the system. The reason is that, as noted above, the 
projector 112 (2/2) is not rank-one, so there is no guarantee 
that it will disentangle the system from the bath. So 
the monitoring procedure Diosi describes cannot possibly 
correspond to the solution of the NMSEE ([7]). Moreover, 
it is easy to verify that, as expected. 



Tri2 



2|1,2 



iyi,y2] 



2|1,2 



{yi,y2] 



dyidy2 ^ P2- 



That is, the measurements described by Diosi disturb 
the evolution of the system so that it no longer obeys the 
original non-Mar kovian dynamics. 

Markovian limit. — There is one case where Diosi's 
monitoring procedure does give a pure-state solution at 
all times which is identical to that which would be ob- 
tained by an all-at-once measurement at that time. This 
is case a — *■ 0, where y„ = x„. That is to say, the in- 
tial bath state is unentangled, and the apparatuses are 
measured locally. In this Markovian limit we find 



*2(2/i,2;2) 



mil 12/2)2 



■02(2/1,2/2! 



(26) 



where the conditional state 



^^2 (2/1, 2/2); is given by 



cexp 



-(^2-2/2)' 



exp 



-(Xi-j/i)2 IV-o). (27) 



This sequence of exponentials can obviously be continued 
indefinitely. The correspondence between the all-at-once 
measurement and Diosi's monitoring here is not surpris- 
ing: in the Markovian limit the interpretation of a SSE in 
terms of continuous monitoring of the bath is well known. 

To conclude, Diosi has introduced an elegant formu- 
lation of non-Markovian evolution using a local (Marko- 
vian) coupling to the bath but an initially non-local (en- 
tangled) bath state. In this formulation, it is simple to 
monitor the bath without affecting the future evolution 
of the system, because each apparatus only interacts with 
the system once. However, to make the conditioned sys- 
tem state pure, it is necessary to measure not only the ap- 
paratuses which have already interacted with the system, 
but also some of those which are yet to interact. Measur- 
ing the latter necessarily introduces noise that will dis- 
turb the future evolution of the system, so that it will not 



reproduce the original non-Markovian evolution on aver- 
age. We show by explicit calculation that the monitoring 
scheme suggested by Diosi does disturb the evolution in 
this manner, and moreover it even fails to produce pure 
conditional system states. 

While it is certainly possible to derive a non-Markovian 
stochastic Schrodinger equation, its solution can only be 
interpreted as a conditioned system state at some par- 
ticular (but arbitrary) time i [y, @| . Connecting the so- 
lutions at different times creates the illusion of a "quan- 
tum trajectory", but is not part of standard quantum 
mechanics. Rather, it is related to Bohmian mechanics 
and its generalizations which also allow one to derive 
discontinuous (jumplike) trajectories Whether the 

jumplike non-Markovian trajectories recently introduced 
in Ref. can be interpreted in a similar manner re- 
mains to be determined. But from the arguments in this 
Letter we know that non-Markovian pure-state trajecto- 
ries cannot be interpreted as true quantum trajectories. 
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